Using error correction to determine the noise model 
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Quantum error correcting codes have been shown to have the ability of making quantum infor- 
mation resilient against noise. Here we show that we can use quantum error correcting codes as 
diagnostics to characterise noise. The experiment is based on a three-bit quantum error correcting 
code carried out on a three-qubit nuclear magnetic resonance (NMR) quantum information pro- 
cessor. Utilizing both engineered and natural noise, the degree of correlations present in the noise 
affecting a two-qubit subsystem was determined. We measured a correlation factor of c = 0.5 ± 0.2 
using the error correction protocol, and c = 0.3 ± 0.2 using a standard NMR technique based on 
coherence pathway selection. Although the error correction method demands precise control, the 
results demonstrate that the required precision is achievable in the liquid-state NMR setting. 

PACS numbers: 03.67.Pp, 05.40.-a,33.25.+k 



I. INTRODUCTION 

The idea of using quantum mechanical systems as in- 
formation processing devices was proposed more than 
two decades ago and yet, experimental realization of 
such devices remains a challenge. Ultimately, all physical 
realizations are faced with the presence of decoherence, 
or noise, caused by uncontrollable interactions with the 
environment Q. 

To prevent the loss of coherence in the quantum me- 
chanical processor, the theory of quantum error correc- 
tion (QEC) has been developed 0,11 US]- QEC 

works by 

encoding the state of a system lying in a certain Hilbcrt 
space into a state in a larger Hilbert space. The encoding 
is designed to make possible the recovery of the original 
information after noise has acted on the overall system, 
through decoding and syndrome measurement, as longas 
the the noise level falls below a certain threshold @, H| ■ 
Many quantum error correction codes (QECC) have been 
developed for specific classes of noise models. For exam- 
ple, there are codes that can correct arbitrary single qubit 
errors 0, 0, H, 0] but fail to correct multi-qubit errors. 
This work shows how this failure can be used to extract 
information about the noise of the system. 

Most QECC are developed for independent, or uncor- 
related, error models, meaning that the errors happening 
on one qubit are assumed to be independent of the errors 
on other qubits. Clearly, knowing whether or not there 
exist correlations in the noise model plays an important 
role in choosing the best QECC for a given system. 

The noise model can be determined exactly by per- 
forming process tomography [n| EH- However, the num- 
ber of experiments required for complete tomography 
grows exponentially in the number of qubits. Often, pro- 
cess tomography is not needed and important (but par- 
tial) information about the noise can be extracted from 
fewer experiments. Here, we demonstrate the use of a 
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three-qubit QECC to extract the two-qubit noise cor- 
relation factor under a transverse relaxation process in 
nuclear magnetic resonance (NMR) (e.g. T2 relaxation). 

Transverse relaxation is the main source of decoherence 
in liquid state NMR. NMR was used to perform the first 
experimental implementation of QEC [12j . where it was 
shown that the three-qubit QECC could correct singlc- 
qubit phase flip errors caused by T2 relaxation. Here, we 
will first briefly review the fundamentals of NMR, then 
model the noise present in such systems and show how 
noise correlations can affect the fidelity of the QEC pro- 
tocol. We will describe a series of natural and engineered 
noise experiments for determining the two-spin noise cor- 
relations present for the 13 C subsystem of acetyl chloride 
(dissolved in deuterated chloroform). The experimental 
results arc in agreement with those obtained using the 
standard NMR technique of coherence selection. In light 
of the exquisite sensitivity of these experiments to con- 
trol imperfections, the results also demonstrate the high 
degree of precision attainable in controlling nuclear spins 
in liquid state NMR. 

The results demonstrate that QEC can not only be 
used for correcting the effects of decoherence, but can also 
help to characterize the nature of those errors. Moreover, 
as QEC is a requirement for scalable quantum informa- 
tion processing (QIP), this methodology is universal for 
probing noise correlations in physical systems suitable for 
QIP. 



II. NMR QUANTUM COMPUTING 

Liquid state nuclear magnetic resonance has proven 
to be a useful system for experimentally benchmark- 
ing small-scale quantum information processing devices 
E2, El, E3, EE El| ■ A NMR quantum information pro- 
cessor consists of an ensemble containing of order 10 20 
molecules with spin-i nuclei dissolved in a liquid sol- 
vent. Placed in a strong homogeneous magnetic field, 
nuclear spins precess about the direction of the field, de- 
fined conventionally as the z-axis. The rate at which the 
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spins precess is the Larmor frequency, and nuclei with 
distinct Larmor (resonance) frequencies can be mapped 
to qubits. In the liquid state, picosecond-scale rotation 
and translation of the molecules causes spins on separate 
molecules to effectively decouple on the NMR timescale. 
Therefore, to a very good approximation, all molecules of 
the same type experience identical environments and the 
Hilbcrt space of the nuclear spin ensemble can be taken 
as that of a single molecule. Moreover, the rotational 
degree of freedom causes the internal dipolar interaction 
between the spins on each molecule to vanish. At thermal 
equilibrium, the Boltzmann distribution gives a slight ex- 
cess of spins pointing along the +z direction, so that an 
average magnetization is present along +z. 

Control of the qubits is achieved by a radio-frequency 
(RF) Hamiltonian in which the frequency, phase and du- 
ration of the RF can be controlled externally. Singlc- 
qubit rotations are performed using RF pulses resonant 
with the Larmor frequency of the targeted qubit. By 
varying the RF duration and phase, rotations of arbitrary 
angle can be generated about any axis in the xy-plane. 
Two-qubit operations additionally use the natural cou- 
pling terms present in the internal Hamiltonian, which 
will be elaborated below. 



III. HAMILTONIAN AND NOISE MODEL 

The nuclear spin Hamiltonian in liquid-state NMR is 
composed of two types of terms, one corresponding to 
the single-spin Zccman interaction (the term leading to 
precession) and bilinear terms corresponding to the scalar 
spin-spin coupling (J-coupling). For a molecule with N 
spin- 1/2 nuclei, the weak coupling Hamiltonian is given 
by 
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1 N 7T 



(1) 



where Vi are the Larmor frequencies, Jy is the coupling 
strength between spins i and j, and Zi is the z Pauli 
matrix for spin i. Note that when the condition — 
Vj\ >> Jij/2 docs not hold (strong coupling regime), 
the scalar coupling operator takes the more general form 
<7i ■ (T j ~ XiXj + Y-iYj + ZiZj. 

Despite the motional averaging that occurs in the liq- 
uid state, the ZZ part of the intermolecular dipolar 
Hamiltonian is capable of creating relaxation, while the 
XX + YY part still averages to zero due to the weak 
coupling approximation 17] . This ZZ interaction couples 
with the molecular motion and give rise to rapidly fluctu- 
ating local magnetic field, which effectively presents itself 
as a variation of the Larmor frequencies. This process is 
known as transversal relaxation. 

Consider a single spin qubit surrounded by an envi- 
ronment £ consisting of N other spins -i. The dipolar 
coupling between the qubit and its environment is de- 



scribed by the unitary evolution 
U 
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where bj is the interaction strength between the qubit 
and the j th spin of the environment for a certain amount 
of time. 

The global system can be assumed to initially be in 
the state 



Pglob = Rini ® RE 



(3) 



where p^ represents the initial state of the qubit and p£ 
is the state of the environment. After the interaction with 
its environment, the final state of the qubit will be given 
by partial tracing the environment system. Moreover, 
the interaction strengths have a certain distribution of 
value q(b), so that the state affected by the noise have 
the form 
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where we have defined £ a = ~}2 m b m {— l) a *™, a m being 
the m th digit of a. In room temperature liquid state 
NMR, the deviation of the state of the environment 
from the completely mixed state is negligible, so that 
(a\p£\a) = l/N. Because the environment is isotropic, 
the distribution q(b) is a symmetric function of the 6j's. 
Therefore, the summation over a in Eq. [3] can be ab- 
sorbed in a new distribution of b and by letting a - 
The final state is then represented by 
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where a = b m and p{a) is the distribution of a which 
takes into account the new distribution of the bj's. The 
interaction of the qubit with the environment causes an 
incoherent averaging of z rotations, which is equivalent to 
a variation of the Larmor frequency of the qubit. In liquid 
state NMR, N is a large number and the central limit 
theorem indicates that a has a gaussian distribution. For 
a M qubit system, this model generalizes to 
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where a — (a±, . . . , oa/), Z = (Z\, . . . Zm) and p(a) i 
the multivariate gaussian distribution [18[ 
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S is the covariant matrix, or the correlation matrix, 
which takes the form 



-'VI — Y^i i 

~2 



(8) 
(9) 
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where of is the variance of cti. Cij is the correlation factor 
between a, and a, , which has value 



(10) 



From the Cauchy-Schwarz inequality, < Cy < 1. For a 
single qubit, such a noise model will affect the state as 

|*X/| -> eT<i-^0| fc )^| (11) 

From empirical results of transverse relaxation in NMR, 
the state of a single spin decays in time as 



,-(l-<5 fc z)72t 



(12) 



where I/72 = ?2 is the relaxation time constant. The 
variance of the distribution of the interaction strength of 
a qubit with its environment can thus be related to its 
relaxation time constant by 



272*. 



(13) 



For two qubits, the noise correlation factor will affect 
the decay of their mutual state as follow: 



\km)(ln\ -> er {1 - &kl) ^ '-(!-«» 
x \ km)(ln\, 
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(14) 



where rjij = |((— l) 1 — (— l)- 7 ). If correlations in the noise 
affecting two qubits is present, the transverse relaxation 
will be faster for a two spin double quantum coherence 
(e.g. |00)(11| and |11)(00|) and slower for a two spin zero 
quantum coherence (e.g. |01)(10| and |10)(01|). 

The correlation in the noise on two qubits can be un- 
derstood through distinguishability. If two nuclei pro- 
cess at the same Larmor frequency, they are magnetically 
equivalent and thus see the same environment. The two 
spins will interact identically with the environment, thus 
yielding a correlation factor of 1. Two spins of different 
nuclear specie are distinguishable and the environment 
will act differently on each of them. No correlation is 
expected the respective noise i.e. e i2 = 0. If we consider 
two nuclei of the same species with slightly different Lar- 
mor frequency, they are distinguishable enough to per- 
form independent control, but they are chemically "near 
indistinguishable" . The effect of the environment is thus 
partially correlated, i.e < en < 1- 



IV. ENGINEERING THE NOISE FOR TWO 
QUBITS 

By explicitly expanding Eq. [7]for two qubits, the noise 
model takes a discrete Kraus form, 



Pi 
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where the unitary Kraus operators £/, and their coeffi- 
cients pi are given in Table |TJ One can thus engineer the 
noise on two qubits with a series of six separate experi- 
ments, each of them implementing a different Kraus op- 
erator, and then adding the results with the correspond- 
ing coefficient. This Kraus decomposition demonstrates 
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TABLE I: Kraus decomposition for the correlated noise on 
two qubits 

that the transversal relaxation in NMR is equivalent to 
a phase flip error, where the qubits undergo a phase flip 
given by the operator in the first column of Table U with 
a probability given by the second column. 



V. DETERMINING THE CORRELATION 
FACTOR 

This section will explain how the noise correlation fac- 
tor between two spins can be extracted using standard 
NMR techniques and how quantum error correction can 
be used to achieve similar results. The details and re- 
sults of the experimental implementation, as well as a 
summary of the advantages of this new technique will 
then conclude this section. 



A. NMR techniques 

In NMR, measurement of transversal relaxation times 
(T2's) is a standard techniques and is implemented 
through single coherence decay and spin echo [Til]. The 
same technique is applicable to double coherences to ex- 
tract the noise correlation factor between two spins. Con- 
sider the following pulse sequence: 



7Tl7r 2 



(16) 



where r is a certain time delay and iii correspond to 
a 7T pulse on nuclei i around any axis in the xy-plane. 
They are used to refocus the field inhomogeneities via 
spin echo. If we apply such a pulse sequence to a state 
of the form 



[ooxiil, 



(17) 
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FIG. 1: The quantum circuit of the three qubit quantum 
error correction code. The two qubit gates are C-NOTs and 
the three qubit gate is a Toffoli gate. 

which can be created using standard NMR techniques 
of coherence selection such as phase cycling [2(| or field 
gradients, the noise model developed earlier predicts that 
the amplitude of such a state should decay as 

|00)(11| -> e ^ 1)t -^ 2) '- 2ci2t V / '^^|00)(ll|. (18) 

In NMR, only single coherence state can be detected. A 
final \ pulse is thus needed on one of the spin to detect 
such a state. By repeating the experiment for various 
value of t, one obtain a decay curve. Once the values 
of Ti are measured using single coherence decay experi- 
ments, it is possible to deduce the value of c\i . 



where it has been assumed that the noise affecting qubit 
3 was uncorrelated with the other qubits (because qubit 
3 is represented by a separate nuclear species in our ex- 
periment). The correlation factor can be extracted from 
the deviation of the fidelity from unity, due to the failure 
of the code. 

As demonstrated in Sec. IIVI it is possible to engineer 
the correlated noise on two spins using six different ex- 
periments. If we want to engineer the noise for a third un- 
correlated qubit, it is done with twelve experiments, us- 
ing the union of two sets of Kraus operators/coefficients 
given by 

{U' k ,p' k } = {Ul' 2 ,(l- q )p]- 2 }\J{U^Z 3 , qP 1 i 2 }(21) 

for k = 1 ... 12 and i, j = 1 ... 6 and where q = |(1 — 
e 72 ) corresponds to the probability of the uncorrelated 
qubit to undergo a phase flip and U i 1 and p\ 2 are the 
correlated noise Kraus operators/coefficients given in Ta- 
ble D 

Therefore, we can implement the QECC using those 
twelve noise operators and obtain the fidelity decay for 
various value of C12 and t . 

C. The experiment 



B. Three qubit quantum error correction code 

The three qubit quantum error correction code can 
protect one qubit of information \ip) against single qubit 
errors about one Pauli axis. The quantum circuit for this 
code can be found in Fig. [TJ If errors happen during the 
noise period, it can be shown that this code corrects any 
single qubit phase error, i.e. errors of the form Z\,Z-i 
or Z%, but fails at correcting multiple phase errors, i.e. 
Z\Z2, Z\Z%, Z2Z3 and ZxZ^Z^. 

As seen above, the natural noise present in NMR con- 
sists of a phase flip. A valid measure to quantify the effect 
of the noise on the system is to consider the fidelity of en- 
tanglement Fe [lH , which corresponds to averaging the 
state-correlation for the density matrix states X, Y and 
Z. In other words, the state-correlation f ^ for an initial 
state A consist on the amount of polarization in the out- 
put relative to the input. The fidelity of entanglement is 
then given by 



1 
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(19) 



If the decoherence is caused solely by the transversal re- 
laxation, the fidelity of entanglement over time of such a 
protocol is given by 



F„ 
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The theory laid down in the previous section assumed 
that the system is composed of two noise correlated 
qubits and one uncorrelated qubit. As seen in section 
III! such a system can be found in a molecule containing 
two spins of the same species with different Larmor fre- 
quency and one of a different kind. For this experiment, 
we have chosen the 13 C-labeled acetyl chloride dissolved 
in deuterated chloroform and used a 700 MHz Brukcr 
Avance NMR spectrometer with dual inverse cryoprobc. 
The structure, chemical shifts and J-coupling strengths 
of the molecule are given in Fig. [5] For this molecule, the 
assumption of weak coupling used throughout section Unl 
is fulfilled due to the large chemical shift difference be- 
tween the two carbons. 

The T^'s for each nuclei have been determined us- 
ing a series of spin echo experiments for various delays 
and their values are given in Table HH To implement 
the quantum error correction code on this molecule, the 
circuit in Fig. [1] was first converted into gates imple- 
mentable in NMR, which consist of single qubit rota- 



Nucleus 


Ti 


T 2 


M 


4.0 ±0.1 


1.2 ±0.1 


Ci 


7.9 ±0.4 


2.1 ±0.1 


Ci 


15.2 ± 0.8 


0.24 ±0.03 



TABLE II: Ti and T2 values for the acetyl chloride. It can 
be seen that for the maximal duration of the experiment (~ 
300 ms), the effect caused by Ti relaxation can be neglected. 
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FIG. 2: The 13 C-labeled acetyl chloride molecule. The diago- 
nal elements of the table gives the chemical shift (difference in 
Larmor frequency) for each nucleus with respect to a reference 
frequency (700.13 MHz for the hydrogens and 176.05 MHz for 
the carbons). The three hydrogens forming the methyl group 
are indistinguishable and form a spin-i and § subspace. The 
spin- 1 subspace was selected using a three-step pulse sequence 
and "crusher" field gradients [22l ]. 



tions about any axis in the xy-plane or around the z- 
axis, and J-coupling evolutions. A J-coupling of length 
t = jj is locally equivalent to a C-NOT. Moreover, the 
z-rotations can be done instantaneously by changing the 
phase of subsequent pulses. This ideal NMR pulse se- 
quence was then fed into a homemade compiler which 
estimates the first order phase and coupling errors dur- 
ing the pulses and then tracks the phase of the subse- 
quent pulses and optimize the refocusing scheme and J- 
coupling delays to minimize overall coupling errors [To] ]. 
Spatial averaging [l5j was used to initialized the states 
|00)(00|i 2 <g> X 3 , |00)(00|i2 <8> Y 3 and |00)(00|i 2 ® Z 3 , from 
the thermal state of a liquid state NMR system. 

The quantum error correction code was first imple- 
mented using engineered errors in twelve experiments. 
The purpose of analyzing engineered noise is to be able 
to generate different fidelity of entanglement decay curves 
corresponding to different value of correlation factors. It 
is done by adding each experiments weighted by the cor- 
responding coefficient given in Table [T] Once a fidelity 
decay curve is obtained for the natural noise, it is pos- 
sible to extract the correlation factor by comparing to 
which engineered noise fidelity decay curve the natural 
noise curve correspond to. 

The natural noise fidelity decay curve was obtained by 
implementing the identity map during the noise section of 
Fig. [T] To perform this implementation the spins could 
not be simply decoupled from one another using multi- 
ple 7T pulses (e.g. the Hadamard refocusing scheme [IH). 
Under such a refocusing scheme for a time t, the double 
coherence terms in a density matrix spend as much time 
in zero quantum coherence as in double quantum coher- 
ence. Therefore, from Eq. 1141 the correlation factor term 
in the exponential decay cancel and do not affect the de- 
cay of the double coherence term. 

If we let the natural noise act on the system for a 



a double coherence for the two carbons have remained 
in double coherence during the entire delay. The field 
inhomogencitics can be rcfocused by applying simultane- 
ous 7T pulses on the carbons which leaves the J-coupling 
evolution untouched. 
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identity and the terms of the density matrix containing 



FIG. 3: Experimental results. The lines correspond to the 
fidelity decay for noise correlation factors of 0, 0.5 and 1 as 
a function of time simulated from the measured T2's and the 
experimental fidelities obtained by implementing engineered 
noise. The points are the fidelities when the system is affected 
by natural noise for a various amount of time. 

The fidelities of our experiments have been ex- 
tracted by fitting every peak of the NMR spectra using 
Lorcntzian shape curves. The resulting values can be 
seen in Fig. [31 where the curves for the fidelity of entan- 
glement for engineered noise are shown for correlation 
factors of 0, 0.5 and 1. Ten values of the fidelity de- 
cay obtained by applying natural noise are shown, from 
which we extracted a correlation factor of C12 = 0.5 ±0.2. 

This experiment needs a high degree of precision, since 
within the time interval used to implement natural noise 
(from to 320 ms) , the maximum difference between the 
C12 = and C12 = 1 curves is 3%. Integrating the square 
of the noise of a spectrum over a region corresponding to 
the width of a signal peak estimates the signal to noise 
ratio to be of the order of 1%. Therefore, the fluctuation 
of the measured fidelities due to the noise explains the 
large uncertainty on the measured correlation factor. 

Using the usual NMR technique of double coherence 
decay, the noise correlation factor between C\ and C2 
was determinded to be c\% = 0.3 ± 0.2. The interval 
using QECC agrees with the value obtained using the 
traditional double-quantum coherence decay technique, 
to within experimental error. 



D. Discussion 

By comparing the above two techniques to extract the 
noise correlation factor between two spins, one could ar- 
gue that the QECC technique is much more involved 
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then the standard NMR technique, while yielding to the 
same conclusion. The goal of the present experiment was 
to demonstrate that the use of QEC to probe the noise 
present in a system is feasible and that the control nec- 
essary to get the error information is achievable. From 
there, it is possible to generalize this technique to any 
physical system with more complex noise model. If the 
noise contains not only phase errors, but also bit errors 
and/or a combination of the two, the same technique 
could be applied using more complex QECC, such as the 
five qubit code @. 

Other technical advantages arise from the signal de- 
tection. Using the NMR technique, there is an doubly 
exponential decay in the signal amplitude for a double 
coherence decay (see Eq. [M]) . From the nature of the 
QECC technique, the signal decays slower, thus allowing 
better statistics and analysis. If the system under analy- 
sis contains three spins of the same type, there would be 
a possibility of three different correlation factors C12, C23 
and C13. Using standard NMR technique, three different 
experiments with different initial states would be needed 
to extract those three values. Using the QECC tech- 
nique, only the noise portion of the pulse sequence would 
need adjustment by changing the rcfocusing scheme in 
order to refocus the unwanted correlation, e.g. a tt pulse 
on qubit 1 would cancel the correlation C12 and C13 for 
the reasons explained earlier. 

Finally, this technique could be used to validate our 
assumption that the noise is effectively gaussian. In the 
case where the system contained three noise correlated 
spins, our gaussian assumption ensures that the noise is 



only pairwise correlated, i.e. C123 = 0. If it is so, a triple 
coherence decay curve should be described using only the 
T2 values and the pairwise correlation factors. This curve 
can be obtain using the NMR technique of triple coher- 
ence decay, but would yield a curve that decays triply 
exponentially. On the other hand, a curve affected by a 
triple correlation factor could be obtained by the same 
QECC pulse sequence by letting all the noise correlations 
act during the noise part of the pulse sequence. In the 
case where that curve would not be described properly 
using only the pairwise correlations, it would be an indi- 
cation of the failure of the noise model. 



VI. CONCLUSION 

In this work, we have demonstrated that QEC can be 
used to probe a physical system and extract partial, but 
important information about the noise model without 
having to perform full quantum process tomography. The 
technique was implemented successfully in a liquid-state 
NMR quantum information processor, but is applicable 
to any QIP device in which standard quantum error cor- 
rection can be carried out. 
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